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Abstract. The Fe isotope effect (Fe-IE) on the transition temperature Tc and the 
crystal structure was studied in the Fe chalcogenide superconductor FeSei_a; by means 
of magnetization and neutron powder diffraction (NPD). The substitution of natural 
Fe (containing ~ 92% of ^^Fe) by its lighter ^"'Fe isotope leads to a shift of Tc of 
0.22(5) K corresponding to an Fe-IE exponent of ape = 0.81(15). Simultaneously, 
a small structural change with isotope substitution is observed by NDP which may 
contribute to the total Fe isotope shift of Tc- 
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Historically, the isotope effect played a crucial role in elucidating the origin of 
the pairing interaction leading to the occurrence of superconductivity. The discovery 
of the isotope effect on the superconducting transition temperature Tc in Hg [1] in 
1950 provided the key experimental evidence for phonon- mediated pairing as formulated 
theoretically by BCS subsequently. The observation of unusually high Tc's in the newly 
discovered Fe-based superconductors immediately raised the question regarding the 
pairing glue and initiated isotope effect studies. Currently, we are aware of two papers 
on isotope experiments with, however, contradicting results. Liu et al. [2] showed that 
in SmFeAsO0.85F0.15 and Bao.6Ko.4Fe2As2 the Fe isotope effect (Fe-IE) exponent, 

«Fe = -dlnT./dlnMpe = -(AT,)/T,)/(AMFe/MFe), (1) 

reaches values of a-pe — 0.35 (Mpe is the Fe atomic mass), while Shirage et al. [3] found a 
negative Fe-IE exponent a-pc — —0.18 in Bai_2;K2.Fe2As2. Note, that the only difference 
between the Bai_a;Ka;Fe2As2 samples studied in Refs. [2] and [3] was the preparation 
procedure (low-pressure synthesis in [2] vs. high-pressure synthesis in [3]), while the 
potassium doping [x ~ 0.4) as well as the Tc's for the samples containing natural Fe 
(Tc ~ 37.3 K in |2] vs. ~ 37.8 K in [3]) were almost the same. 

In this paper we study the Fe-IE on and on the structural parameters (such as 
the lattice parameters a, b, and c, the lattice volume V, and the distance between the 
Se atom and Fe plane, Se height /ise) for another representative of the Fe-based high- 
temperature superconductors (HTS), namely FeSei_a;. The substitution of natural Fe 
(containing ~ 92% of ^^Fe) by its lighter ^^Fe isotope leads to a shift of Tc of 0.22(5) K 
corresponding to an Fe-IE exponent of ctpe = 0.81(15). 

The ^"^FeSei-aj/^^FeSei-a; samples (here after we denote natural Fe containing 
~ 92% of ^^Fe isotope as ^^Fe) with the nominal composition FeSeo.gs were prepared by 
a solid state reaction made in two steps. Pieces of Fe (natural Fe: 99.97% minimum 
purity, average atomic mass Mpe = 55.85 g/mol, or ^"^Fe: 99.99% purity, 99.84% isotope 
enriched, Msipf. = 54.0 g/mol) and Se (99.999% purity) were first sealed in double walled 
quartz ampules, heated up to 1075°C, annealed for 72 h at this temperature and 48 h at 
420°C, and then cooled down to room temperature at a rate of 100°C/h. As a next step, 
the samples, taken out of the ampules, were powderised, pressed into pellets, sealed into 
new ampules and annealed first at 700°C for 48 h and then at 400°C for 36 h, followed 
by cooling to room temperature at a rate of 200°C/h. Due to the extreme sensitivity 
of FeSei_a; to oxygen [1], all the intermediate steps (grinding and pelletizing) as well 
as the preparation of the samples for the neutron powder diffraction and magnetization 
experiments were performed in a glove box under He atmosphere. 

The Fe-IE on the structural properties was studied by neutron powder diffraction 
(NPD) experiments by using the high-resolution powder diffractometer HRPT (Paul 
Scherrer Institute, Switzerland) [5]. The experiments were carried out at a wavelength 
A = 1.494 A. The ^^FeSei-^^/^^FeSei-i: samples, placed into vanadium containers, were 
mounted into a He-4 cryostat in order to reach temperatures between 5 and 250 K. High 
statistics data were taken at 250 and 5 K. Data at 10 < T < 240 K were collected with 
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Figure 1. (Color online) The Rietveld refinement pattern and difference plot of NPD 
data for ^^FeSei.^; (panel a) '^^FeSei.j; (panel b) at T = 250 K. The rows of ticks show 
the Bragg-peak positions for the main phase FeSe (PAnmm) and two impurity phases: 
Fe {ImSm) and hexagonal FeSe {PQ^/mmc). The main tetragonal phase corresponds 
to 0.975(5) and 0.975(4) Se occupancy for ^^FeSei-^. and ^""FeSei-a;, respectively. 

Figure [T] shows the NPD spectra taken at T = 250 K. The differences in peak 
intensities, clearly visible at small 6, are caused by the different values of the coherent 
neutron scattering length (bcoh) of natural Fe and that of the ^'^Fe isotope. The 
refinement of the crystal structure was performed by using the FULLPROF program 
[6] with b^^l = 9.45 • 10"^^ m, bj^^" = 4.2 • 10"^^ m, and b^^j, = 7.97 ■ 10"^^ m [7]. The 
refined structural parameters at T = 250 K and 5 K are summarized in Table 1. The 
amount of the impurity phases and the Se content (1 — x), determined for the data sets 
taken at T = 250 K, were kept fixed during the refinement of the NPD spectra at lower 
temperatures. The mass fractions of impurity phases, the hexagonal FeSe {PG^/mmc) 
and Fe (/m3m), were found to be 0.50(10)%, 0.31(4)% and 1.13(18)%, 1.06(7)% for 
^^FeSei_a; and ^^FeSei-^,, respectively. 

Figure [2] shows the temperature dependence of the lattice parameters a, b, and c, the 
lattice volume V, and the Se height hsc of a representative ^''FeSei-a; and a representative 
^^FeSei_a; sample (see Fig. [3]). From Fig. [2^ it is obvious that at ^ 100 K a transition 
from a tetragonal to an orthorhombic structure takes place, analogous to that reported 
in [U |8]. The Fe-IE on the structural transition temperature Tg could be estimated 
from the shift of the interception point of the linear fits to a(T) and b{T) in the 
vicinity of T^, as denoted by the arrows in the inset of Fig. [2^, which was found to 
be AT, = 0.2(2.5) K. Within the whole temperature range (5 K< T <250 K) the 
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Table 1. Structural parameters of ^"^FeSei-^ and ^"^FeSei-^; at T = 250 and 5 K. 
Space group PA/nmm (no. 129), origin choice 2: Fe in (26) position (1/4, 3/4, 1/2); 
Se in (2c) position (1/4, 1/4, z). Space group Cmma (no. 67): Fe in (45) position 
(1/4, 0, 1/2), Se in [Ag) position (0, 3/4, z). The atomic displacement parameters [B) 
for Fe and Se were constrained to be the same. The Bragg R factor is given for the 
main phase; the other reliability factors are given for the whole refinement. 
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lattice constants a and h are slightly larger for ^^FeSei_a: than those for ^^FeSei-^; (see 
Fig. [2^). This is in contrast to the lattice parameter c, which within the same range 
is marginally smaller for ^^FeSei_j; than for ^^FeSei_a; (Fig. Wp)- The lattice volume 
remains, however, unchanged. Consequently, substitution of ^^Fe by ^^Fe leads to a 
small, but detectable enhancement of the lattice along the crystallographic a and h 
directions and a compression of it along the c— axis, resulting in a change of the shape 
of the Fe4Se pyramid, which is known to influence T^. in Fe-based HTS [H [TOl [11]. This 
is shown in Fig. |2fc where below 100 K the Se atom is located closer to the Fe plane in 
^^FeSei_a; than in ^^FeSei-^,. The corresponding change of the Fe4Se pyramid is shown 
schematically in the inset of Fig. It is important to note that the observed Fe-IE's 
on the lattice parameters are intrinsic and not just a consequence of slightly different 
samples. As shown in Ref. |4], various samples of ^^FeSei_a; with 1 — x ^ 0.98 and 
Tc ^ 8.2 K indeed exhibit the same lattice parameters within experimental error. 

The Fe-IE on the transition temperature T^. was studied by means of magnetization 
experiments. Measurements were performed by using a SQUID magnetometer 
(Quantum Design MPMS-7) in a field of ixqH = 0.1 mT for temperatures ranging from 
2 to 20 K. In order to avoid artifacts and systematic errors in the determination of the 
isotope shift of Tc it is important to perform a statistical study: i.e. to investigate 
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series of ^^FeSei_a;/^^FeSei_a; samples synthesized exactly the same way (the same 
thermal history, the same amount of Se in the initial composition). The magnetization 
experiments were conducted for six ^^FeSei-^: and seven ^^FeSei_a; samples, respectively. 
The inset in Fig. [3] shows an example of zero-field cooled (ZFC) magnetization curves 
for a pair of ^^FeSei-^/^^FeSei-i; samples {Mnorm was obtained after subtracting the 
small paramagnetic offset M^agn measured at T > Tc and further normalization of 
the obtained curve to the value at T = 2 K, see Fig. 1 in Ref. [1] for details). The 
magnetization curve for ^^FeSei-^,. is shifted almost parallel to higher temperature, 
implying that of ^^FeSei.j; is higher than that of ^^FeSei-^;. The resulting transition 
temperatures determined from the intercept of the linearly extrapolated Mnorm{T) 
curves with the M = line for all samples investigated are summarized in Fig. |3l 
The Tc's for both sets of ^^FeSei_a;/^^FeSei_a; samples fall into two distinct regions: 
8.39 < ^^Tc < 8.48 K and 8.15 < ^^Tc < 8.31 K, respectively. The corresponding mean 
values are: ^^T^ = 8.43(3) K and ^^T^ = 8.21(4) K. Note, that one out of the seven 
^^FeSei_a; samples had ~ 8.44 K which is by more than 5 standard deviations above 
the average calculated for the rest of the six samples. We have no explanation for this 
discrepancy, but decided to show this point for completeness of the data collected. 

The Fe-IE exponent apc was determined from the data presented in Fig. [3] using 
Eq. dl]), where the relative Fe isotope shift of the quantity X is defined as AX/X = 
{^^X X) (this definition of AX/X is used throughout the paper). With 
s^Tc = 8.43(3) K, s^T, = 8.21(4) K, M54fc = 54 g/mol, and Msepe = 55.85 g/mol 
one obtains apc = 0.81(15). Two points should be emphasized: i) The positive sign of 
the Fe-IE exponent ape is similar to that observed in phonon mediated superconductors, 
such as elemental metals P and MgB2 [12] as well as in cuprate HTS [131 E] where the 
pairing mechanism is still under debate. Bearing in mind that a positive Fe-IE exponent 
was also observed in SmFeAsO0.85F0.15 and Bao.6Ko.4Fe2As2 [2], we may conclude that 
at least for three compounds representing different families of Fe-based HTS (1111, 122, 
and 11) the sign of the Fe-IE on is conventional. This suggests that the lattice plays 
an essential role in the pairing mechanism in the Fe-based HTS. ii) The Fe-IE exponent 
a-pe = 0.81(15) is larger than the BCS value a^'^^ = 0.5 as well as more than twice 
as large as ape — 0.35 reported for SmFeAsO0.85F0.15 and Bao.6Ko.4Fe2As2 [2]. Note 
that an enhanced value of the oxygen isotope exponent (ao — 1) was also observed 
in underdoped cuprate HTS [H] and was shown to be a consequence of the polaronic 
nature of the supercarriers in that class of materials [I5]. Recently, Bussmann- Holder 
et al. [16] showed that in the framework of a two-band model polaronic coupling in the 
larger gap channel as well as in the interband interaction induce a (doping) dependent 
Fe-IE: ape increases strongly with reduced T^. (doping), reaching ctpe — 0.9 at Tc ^ 10 K. 
Note that a similar generic trend is observed in cuprate HTS [T3l [T3] . 

However, our structural refined NPD data suggest that part of the large Fe-IE 
ape = 0.81(15) may result from the tiny structural changes due to ^^Fe/^^Fe substitution. 
In the following we discuss a possible structural effect on the observed Fe-IE on Tc. It 
is known that in FeSei_3; a decrease of the Se height caused by compression of the 
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Fe4Se pyramid leads to an increase of T, by AT^^^- /{Ahsc/hsc) ^ 3.4 K/% dUIII]. In 
contrast, an increase of the Se(Te)-Fe-Se(Te) angle in the FeSei_j^Tej^ family (angle (3 
in our notation [18], see the inset of Fig. [2t) results for y < 0.5 in a decrease of Tc by 
ATf/(A/3//3) ~ 2.9 K/% [10]. Based on the structural data presented in Fig. [2] one 
obtains Ahse/hsc = 0.22(8)% and A(3/(3 = -0.13(4)%, leading to AT> = 0.7(3) K 
and AT^ = —0.4(2) K (in this estimate the values of the lattice constants a and b, 
and hse were averaged over the temperature regions denoted as solid lines in Figs. [2^ 
and c). Bearing in mind that all Fe-based HTS are similarly sensitive to structural 
changes as FeSei-^^ (see, e.g., [9l [TOl [TTl [19]) we conclude that the shift of caused 
by tiny modifications of the crystal structure upon isotope exchange may contribute to 
the total Fe-IE exponent. However, the large errors of AT^^'= and ATj^ do not allow a 
reliable estimate of this structural effect on the Fe-IE on Tc. A better estimate of this 
effect can be made based on the empirical relation between Tc and the lattice parameter 
a for the 11 family FeSci-yTe^ [TUl US]. Assuming that the relation Tc(a) is also valid 
for FeSei_a; one obtains from the data presented in Ref. [I9] for y < 0.5 the relation 
AT^/{Aa/a) ^ 6 K/%. With (Aa + Ab)/{a + b) = 0.0195(14)% this gives rise to a 
structural shift of Tc of AT^^^ ^ 0.1 K (the lattice constants a and b were averaged over 
the temperature regions marked as solid line in Fig. [2^). Taking this correction into 
account yields a rough estimate of the intrinsic Fe-IE exponent of a^^ ~ 0.4. This value 
is comparable with apc — 0.35 reported for SmFeAsO0.85F0.15 and Bao.6Ko.4Fe2As2 [2]- 

To summarize, the ^^Fe/^^Fe isotope effects on the superconducting transition 
temperature and the crystal structure were studied in the iron chalcogenide 
superconductor FeSei.^;. The following results were obtained: (i) The substitution 
of the natural Fe (Mpe = 55.85 g/mol) by the ^^Fe isotope (M54fc = 54.0 g/mol) gives 
rise to a pronounced Fe isotope shift of the transitions temperature Tc as determined 
by magnetization measurements. The average Tc is found to be ~ 0.22 K higher for the 
^"^FeSei-a; samples as compared to the ^^FeSei_a; samples resulting in a Fe-IE exponent 
of apc = 0.81(15). (ii) The ^^Fe/^^Fe isotope substitution leads to an enhancement 
of the lattice constants a and b and a shrinkage of the lattice constant c. These 
modifications do not affect the lattice volume, (iii) The tetragonal to orthorhombic 
structural transition temperature (T, ~ 100 K) is the same for both ^^FeSei_a; and 
^^FeSei_a; within the accuracy of the experiment, (iv) For temperatures below 100 K 
the distance between the Se atom and Fe plane (Se height) is smaller for the samples 
with ^^Fe. This, together with the results of point ii), imply that ^^Fe/^^Fe isotope 
substitution leads to a compression of the Fe4Se pyramid along the crystallographic 
c— axis and an enhancement along the a— and fe— directions, (v) The structural changes 
caused by ^^Fe/^^Fe isotope substitution induce a shift in Tc which may reduce the 
value of Fe-IE exponent to ^ 0.4, in fair agreement with ape — 0.35 obtained for 
SmFeAsO0.85F0.15 and Bao.6Ko.4Fe2As2 [2]. 

In conclusion, from magnetization experiments the Fe-IE exponent of Tc for the 
FeSei_a: system was determined to be ape = 0.81(15). The tiny changes of the structural 
parameters caused by isotope substitution may contribute to the total Fe-IE exponent. 
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and may help to clarify or even be the origin of the previously reported controversial 
results [21 [3]. However, more detailed and systematic structural investigations on Fe 
isotope substituted samples are required in order to draw definite conclusions. Our 
findings, on the other hand, clearly show that a conventional isotope effect on Tc is 
present which highlights the role of the lattice in the pairing mechanism in this new 
material class. 
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Figure 2. (Color online) The temperature dependence of the lattice constants a and 
h (panel a, in the tetragonal phase a is multiplied by v^), lattice constant c (panel 
b), and the distance between the Se atom and Fe plane /iso (panel c) for ^'^FeSep. 975(5) 
and ^^FeSco. 975(4) samples. The inset in panel a shows the extended part of a{T) 
and b{T) in the vicinity of the structural transition temperature T,, together with the 
linear fits. The inset in panel b represents the temperature dependence of the lattice 
volume V . The inset in panel c shows schematically the modification of the Fe4Se 
pyramid caused by ^^Fe/^'''Fe isotope substitution. The arrows show the direction of 
atom displacements (see text for details). 
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Figure 3. (Color online) The superconducting transition temperature Tc as a function 
of Fe atomic mass for ^^FeSei_a;/^^FeSei_a; samples studied in the present work. The 
open symbols correspond to the samples studied by NPD experiments. The T^s fall 
into the regions marked by the colored stripes with the corresponding mean values 
^^Tc = 8.43(3) K and ^^T^ = 8.21(4) K. The inset shows the normalized ZFC 
magnetization curves Mnorm{T) for one pair of ^^FeSei_a; and ^^FeSei-^^ samples. The 
transition temperature is determined as the intersection of the linearly extrapolated 
MnomriT) curve in the vicinity of Tc with the M — line. 



